Analysis of the limited M. tuberculosis accessory genome reveals potential pitfalls of pan-genome analysis approaches.

Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety of methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. To quantify sources of bias and error related to common pan-genome analysis approaches, we evaluated different approaches applied to curated collection of 151 Mycobacterium tuberculosis ( Mtb ) isolates. Mtb is characterized by its clonal evolution, absence of horizontal gene transfer, and limited accessory genome, making it an ideal test case for this study. Using a state-of-the-art graph-genome approach, we found that a majority of the structural variation observed in Mtb originates from rearrangement, deletion, and duplication of redundant nucleotide sequences. In contrast, we found that pan-genome analyses that focus on comparison of coding sequences (at the amino acid level) can yield surprisingly variable results, driven by differences in assembly quality and the softwares used. Upon closer inspection, we found that coding sequence annotation discrepancies were a major contributor to inflated Mtb accessory genome estimates. To address this, we developed panqc, a software that detects annotation discrepancies and collapses nucleotide redundancy in pan-genome estimates. When applied to Mtb and E. coli pan-genomes, panqc exposed distinct biases influenced by the genomic diversity of the population studied. Our findings underscore the need for careful methodological selection and quality control to accurately map the evolutionary dynamics of a bacterial species.

Self-Hybridization in Leishmania major.

Genetic exchange between different Leishmania strains in the sand fly vector has been experimentally demonstrated and is supported by population genetic studies. In nature, opportunities for Leishmania interstrain mating are restricted to flies biting multiply infected hosts or through multiple bites of different hosts. In contrast, self-mating could occur in any infected sand fly. By crossing two recombinant lines derived from the same Leishmania major strain, each expressing a different drug-resistance marker, self-hybridization in L. major was confirmed in a natural sand fly vector, Phlebotomus duboscqi, and in frequencies comparable to interstrain crosses. We provide the first high resolution, whole-genome sequencing analysis of large numbers of selfing progeny, their parents, and parental subclones. Genetic exchange consistent with classical meiosis is supported by the biallelic inheritance of the rare homozygous single nucleotide polymorphisms (SNPs) that arose by mutation during the generation of the parental clones. In contrast, heterozygous SNPs largely failed to be transmitted in Mendelian ratios for reasons not understood. SNPs that were heterozygous in both parents, however, recombined to produce homozygous alleles in some hybrids. For trisomic chromosomes present in both parents, transmittal to the progeny was only altered by self-hybridization, involving a gain or loss of somy in frequencies predicted by a meiotic process. Whole-genome polyploidization was also observed in the selfing progeny. Thus, self-hybridization in Leishmania, with its potential to occur in any infected sand fly, may be an important source of karyotype variation, loss of heterozygosity, and functional diversity. IMPORTANCE Leishmania are parasitic protozoa that cause a wide spectrum of diseases collectively known as the leishmaniases. Sexual reproduction in Leishmania has been proposed as an important source of genetic diversity and has been formally demonstrated to occur inside the sand fly vector midgut. Nevertheless, in the wild, opportunities for genetic exchange between different Leishmania species or strains are restricted by the capacity of different Leishmania strains to colonize the same sand fly. In this work, we report the first high resolution, whole-genome sequence analysis of intraclonal genetic exchange as a type of self-mating in Leishmania. Our data reveal that self-hybridization can occur with comparable frequency as interstrain mating under experimental lab conditions, leading to important genomic alterations that can potentially take place within every naturally infected sand fly.

Patterns of genomic interrelatedness of publicly available samples in the TB portals database.

The TB Portals program is an international collaboration for the collection and dissemination of tuberculosis data from patient cases focused on drug resistance. The central database is a patient-oriented resource containing both patient and pathogen clinical and genomic information. Herein we provide a summary of the pathogen genomic data available through the TB Portals and show one potential application by examining patterns of genomic pairwise distances. Distributions of pairwise distances highlight overall patterns of genome variability within and between Mycobacterium tuberculosis phylogenomic lineages. Closely related isolates (based on whole-genome pairwise distances and time between sample collection dates) from different countries were identified as potential evidence of international transmission of drug-resistant tuberculosis. These high-level views of genomic relatedness provide information that can stimulate hypotheses for further and more detailed research.

Prokaryotic Genome Annotation.

In the last decade, the high-throughput and relatively low cost of short-read sequencing technologies have revolutionized prokaryotic genomics. This has led to an exponential increase in the number of bacterial and archaeal genome sequences available, as well as corresponding increase of genome assembly and annotation tools developed. Together, these hardware and software technologies have given scientists unprecedented options to study their chosen microbial systems without the need for large teams of bioinformaticists or supercomputing facilities. While these analysis tools largely fall into only a few categories, each may have different requirements, caveats and file formats, and some may be rarely updated or even abandoned. And so, despite the apparent ease in sequencing and analyzing a prokaryotic genome, it is no wonder that the budding genomicist may quickly find oneself overwhelmed. Here, we aim to provide the reader with an overview of genome annotation and its most important considerations, as well as an easy-to-follow protocol to get started with annotating a prokaryotic genome.

Genetic mechanisms of Coxiella burnetii lipopolysaccharide phase variation.

Coxiella burnetii is an intracellular pathogen that causes human Q fever, a disease that normally presents as a severe flu-like illness. Due to high infectivity and disease severity, the pathogen is considered a risk group 3 organism. Full-length lipopolysaccharide (LPS) is required for full virulence and disease by C. burnetii and is the only virulence factor currently defined by infection of an immunocompetent animal. Transition of virulent phase I bacteria with smooth LPS, to avirulent phase II bacteria with rough LPS, occurs during in vitro passage. Semi-rough intermediate forms are also observed. Here, the genetic basis of LPS phase conversion was investigated to obtain a more complete understanding of C. burnetii pathogenesis. Whole genome sequencing of strains producing intermediate and/or phase II LPS identified several common mutations in predicted LPS biosynthesis genes. After passage in broth culture for 30 weeks, phase I strains from different genomic groups exhibited similar phase transition kinetics and elevation of mutations in LPS biosynthesis genes. Targeted mutagenesis and genetic complementation using a new C. burnetii nutritional selection system based on lysine auxotrophy confirmed that six of the mutated genes were necessary for production of phase I LPS. Disruption of two of these genes in a C. burnetii phase I strain resulted in production of phase II LPS, suggesting inhibition of the encoded enzymes could represent a new therapeutic strategy for treatment of Q fever. Additionally, targeted mutagenesis of genes encoding LPS biosynthesis enzymes can now be used to construct new phase II strains from different genomic groups for use in pathogen-host studies at a risk group 2 level.

Draft Genome Sequences of the Avirulent Coxiella burnetii Dugway 7D77-80 and Dugway 7E65-68 Strains Isolated from Rodents in Dugway, Utah.

Here, we present the draft genome sequences of the Coxiella burnetii Dugway 7D77-80 and Dugway 7E65-68 strains, which were isolated from rodents in Dugway, UT, in the 1950s. The strains reside in a distinct genomic group of C. burnetii and are considered avirulent despite having the largest genomes of the Coxiella genus.

Draft Genome Sequences of Historical Strains of Coxiella burnetii Isolated from Cow's Milk and a Goat Placenta.

Here, we report draft genome sequences of historical strains of Coxiella burnetii derived from cow's milk and the placenta of a goat that had aborted. The California and Ohio milk strains display a different sequence type than do contemporary milk strains.

Draft Genome Sequences of Three Coxiella burnetii Strains Isolated from Q Fever Patients.

In the current study, we determined the draft genome sequences of three Coxiella burnetii human disease isolates. The Coxiella burnetii Turkey (RSA315) and Dyer (RSA345) strains were isolated from acute Q fever patients, while the Ko (Q229) strain was isolated from a Q fever endocarditis patient.

Genomic and phenotypic characterization of in vitro-generated Chlamydia trachomatis recombinants.

BACKGROUND: Pre-genomic and post-genomic studies demonstrate that chlamydiae actively recombine in vitro and in vivo, although the molecular and cellular biology of this process is not well understood. In this study, we determined the genome sequence of twelve Chlamydia trachomatis recombinants that were generated in vitro under antibiotic selection. These strains were used to explore the process of recombination in Chlamydia spp., including analysis of candidate recombination hotspots, and to correlate known C. trachomatis in vitro phenotypes with parental phenotypes and genotypes. RESULTS: Each of the 190 examined recombination events was the product of homologous recombination, and no candidate targeting motifs were identified at recombination sites. There was a single deletion event in one recombinant progeny that resulted in the removal of 17.1 kilobases between two rRNA operons. There was no evidence for preference for any specific region of the chromosome for recombination, and analyses of a total of over 200 individual recombination events do not provide any support for recombination hotspots in vitro. Two measurable phenotypes were analyzed in these studies. First, the efficiency of attachment to host cells in the absence of centrifugation was examined, and this property segregated to regions of the chromosome that carry the polymorphic membrane protein (Pmp) genes. Second, the formation of secondary inclusions within cells varied among recombinant progeny, but this did not cleanly segregate to specific regions of the chromosome. CONCLUSIONS: These experiments examined the process of recombination in C. trachomatis and identified tools that can be used to associate phenotype with genotype in recombinant progeny. There were no data supporting the hypothesis that particular nucleotide sequences are preferentially used for recombination in vitro. Selected phenotypes can be segregated by analysis of recombination, and this technology may be useful in preliminary analysis of the relationship of genetic variation to phenotypic variation in the chlamydiae.

Analysis of the genome of leporid herpesvirus 4.

The genome of a herpesvirus highly pathogenic to rabbits, leporid herpesvirus 4 (LHV-4), was analyzed using high-throughput DNA sequencing technology and primer walking. The assembled DNA sequences were further verified by restriction endonuclease digestion and Southern blot analyses. The total length of the LHV-4 genome was determined to be about 124 kb. Genes encoded in the LHV-4 genome are most closely related to herpesvirus of the Simplexvirus genus, including human herpesviruses (HHV-1 and HHV-2), monkey herpesviruses including cercopithicine (CeHV-2 and CeHV-16), macacine (McHV-1), bovine herpesvirus 2 (BHV-2), and a lineage of wallaby (macropodid) herpesviruses (MaHV-1 and -2). Similar to other simplexvirus genomes, LHV-4 has a high overall G+C content of 65-70% in the unique regions and 75-77% in the inverted repeat regions. Orthologs of ICP34.5 and US5 were not identified in the LHV-4 genome. This study shows that LHV-4 has the smallest simplexvirus genome characterized to date.

Footprint of positive selection in Treponema pallidum subsp. pallidum genome sequences suggests adaptive microevolution of the syphilis pathogen.

In the rabbit model of syphilis, infection phenotypes associated with the Nichols and Chicago strains of Treponema pallidum (T. pallidum), though similar, are not identical. Between these strains, significant differences are found in expression of, and antibody responses to some candidate virulence factors, suggesting the existence of functional genetic differences between isolates. The Chicago strain genome was therefore sequenced and compared to the Nichols genome, available since 1998. Initial comparative analysis suggested the presence of 44 single nucleotide polymorphisms (SNPs), 103 small (≤3 nucleotides) indels, and 1 large (1204 bp) insertion in the Chicago genome with respect to the Nichols genome. To confirm the above findings, Sanger sequencing was performed on most loci carrying differences using DNA from Chicago and the Nichols strain used in the original T. pallidum genome project. A majority of the previously identified differences were found to be due to errors in the published Nichols genome, while the accuracy of the Chicago genome was confirmed. However, 20 SNPs were confirmed between the two genomes, and 16 (80.0%) were found in coding regions, with all being of non-synonymous nature, strongly indicating action of positive selection. Sequencing of 16 genomic loci harboring SNPs in 12 additional T. pallidum strains, (SS14, Bal 3, Bal 7, Bal 9, Sea 81-3, Sea 81-8, Sea 86-1, Sea 87-1, Mexico A, UW231B, UW236B, and UW249C), was used to identify "Chicago-" or "Nichols -specific" differences. All but one of the 16 SNPs were "Nichols-specific", with Chicago having identical sequences at these positions to almost all of the additional strains examined. These mutations could reflect differential adaptation of the Nichols strain to the rabbit host or pathoadaptive mutations acquired during human infection. Our findings indicate that SNPs among T. pallidum strains emerge under positive selection and, therefore, are likely to be functional in nature.

Resistance to a novel antichlamydial compound is mediated through mutations in Chlamydia trachomatis secY.

A novel and quantitative high-throughput screening approach was explored as a tool for the identification of novel compounds that inhibit chlamydial growth in mammalian cells. The assay is based on accumulation of a fluorescent marker by intracellular chlamydiae. Its utility was demonstrated by screening 42,000 chemically defined compounds against Chlamydia caviae GPIC. This analysis led to the identification of 40 primary-hit compounds. Five of these compounds were nontoxic to host cells and had similar activities against both C. caviae GPIC and Chlamydia trachomatis. The inhibitory activity of one of the compounds, (3-methoxyphenyl)-(4,4,7-trimethyl-4,5-dihydro-1H-[1,2]dithiolo[3,4-C]quinolin-1-ylidene)amine (MDQA), was chlamydia specific and was selected for further study. Selection for resistance to MDQA led to the generation of three independent resistant clones of C. trachomatis. Amino acid changes in SecY, a protein involved in Sec-dependent secretion in Gram-negative bacteria, were associated with the resistance phenotype. The amino acids changed in each of the resistant mutants are located in the predicted central channel of a SecY crystal structure, based on the known structure of Thermus thermophilus SecY. These experiments model a process that can be used for the discovery of antichlamydial, anti-intracellular, or antibacterial compounds and has led to the identification of compounds that may have utility in both antibiotic discovery and furthering our understanding of chlamydial biology.

Complete genome sequence and annotation of the Treponema pallidum subsp. pallidum Chicago strain.

In syphilis research, the Nichols strain of Treponema pallidum, isolated in 1912, has been the most widely studied. Recently, important differences among T. pallidum strains emerged; therefore, we sequenced and annotated the Chicago strain genome to facilitate and encourage the use of this strain in studying the pathogenesis of syphilis.

Genome sequencing of recent clinical Chlamydia trachomatis strains identifies loci associated with tissue tropism and regions of apparent recombination.

The human pathogen Chlamydia trachomatis exists as multiple serovariants that have distinct organotropisms for different tissue sites. Culture and epidemiologic data have demonstrated that serovar G is more prevalent, while serovar E is less prevalent, for rectal isolates from men having sex with men (MSM). The relative prevalence of these serovars is the opposite for isolates from female cervical infections. In contrast, the prevalence of serovar J isolates is approximately the same at the different tissue sites, and these isolates are the only C-class strains that are routinely cultured from MSM populations. These correlations led us to hypothesize that polymorphisms in open reading frame (ORF) sequences correlate with the different tissue tropisms of these serovars. To explore this possibility, we sequenced and compared the genomes of clinical anorectal and cervical isolates belonging to serovars E, G, and J and compared these genomes with each other, as well as with a set of previously sequenced genomes. We then used PCR- and restriction digestion-based genotyping assays performed with a large collection of recent clinical isolates to show that polymorphisms in ORFs CT144, CT154, and CT326 were highly associated with rectal tropism in serovar G isolates and that polymorphisms in CT869 and CT870 were associated with tissue tropism across all serovars tested. The genome sequences collected were also used to identify regions of likely recombination in recent clinical strains. This work demonstrated that whole-genome sequencing along with comparative genomics is an effective approach for discovering variable loci in Chlamydia spp. that are associated with clinical presentation.

Identification of concomitant infection with Chlamydia trachomatis IncA-negative mutant and wild-type strains by genomic, transcriptional, and biological characterizations.

Clinical isolates of Chlamydia trachomatis that lack IncA on their inclusion membrane form nonfusogenic inclusions and have been associated with milder, subclinical infections in patients. The molecular events associated with the generation of IncA-negative strains and their roles in chlamydial sexually transmitted infections are not clear. We explored the biology of the IncA-negative strains by analyzing their genomic structure, transcription, and growth characteristics in vitro and in vivo in comparison with IncA-positive C. trachomatis strains. Three clinical samples were identified that contained a mixture of IncA-positive and -negative same-serovar C. trachomatis populations, and two more such pairs were found in serial isolates from persistently infected individuals. Genomic sequence analysis of individual strains from each of two serovar-matched pairs showed that these pairs were very similar genetically. In contrast, the genome sequence of an unmatched IncA-negative strain contained over 5,000 nucleotide polymorphisms relative to the genome sequence of a serovar-matched but otherwise unlinked strain. Transcriptional analysis, in vitro culture kinetics, and animal modeling demonstrated that IncA-negative strains isolated in the presence of a serovar-matched wild-type strain are phenotypically more similar to the wild-type strain than are IncA-negative strains isolated in the absence of a serovar-matched wild-type strain. These studies support a model suggesting that a change from an IncA-positive strain to the previously described IncA-negative phenotype may involve multiple steps, the first of which involves a translational inactivation of incA, associated with subsequent unidentified steps that lead to the observed decrease in transcript level, differences in growth rate, and differences in mouse infectivity.